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New results on normalised inclusive jet, dijet and trijet differential cross sections in neutral 
current DIS based on a regularised unfolding procedure are presented. Compared to a 
previously published result on normalised multi-jet cross sections, the new features are an 
extended range in jet pseudorapidity, an improved hadronic energy scale uncertainty of 
1% and the application of an unfolding procedure. The normalised jet cross sections are 
compared to QCD calculations at NLO. The value of the strong coupling determined from 
the normalised inclusive jet, dijet and trijet measurements simultaneously is as{Mz) = 
0.1163 ± 0.0008 (exp.) ±0.0011 (sys.)lg:g2^^ (theo.) ±0.0014 (PDF). 



1 Introduction 

Jet production in neutral current DIS at HERA provides an important testing ground for QCD. 
The measurement of jet quantities is directly sensitive to the strong coupling as and can give 
constraints on the gluon density in the proton. Furthermore, it is a valuable benchmark process 
for Monte Carlo event generators, particularly with regard to parton showers. Two different 
kinds of jet measurements can be distinguished. Inclusive jet measurements, where each single 
jet is counted, and jet measurements like dijet and trijet measurements, where each event that 
fulfills topological and kinematic criteria on jet quantities contributes to the cross section once. 
Both approaches allow to extract the strong coupling by comparing to perturbative QCD pre- 
dictions. 

The measurement presented here is based on data with an integrated luminosity of 361 pb~^ 
collected in the years 2003 - 2007 with the HI detector |lj . The data are identical to a previous 
HI analysis of absolute jet cross sections [2]. This analysis is extended to normalised cross sec- 
tions, where the normalisation is performed with respect to the NC DIS cross section. In ^ and 
in the analysis reported here, improvements on the reconstruction of tracks and the calorimetric 
energy were applied [5]. The correction of detector effects to determine the particle level cross 
section is performed using a regularised unfolding procedure and is presented in more detail in 
this document. 
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2 Phase space definition 



The NC DIS events are selected by requiring an identified scattered electron, a virtuality of 
the exchanged boson (j/TP) of 150 < < 15000 GeV^ and an inelasticity of the interac- 
tion of 0.2 < y < 0.7. The jet finding is performed in the Breit frame of reference where the 
exchanged boson is completely space-like. Particle candidates of the hadronic final state are 
clustered into jets using the inclusive fcx algorithm with a distance parameter Rq = 1, as 
implemented in the Fast Jet package 0. The jets are required to be in the pseudorapidity range 
in the laboratory rest frame between —1.0 < rf^^^^ < 2.5, and the jet momentum in the Breit 
frame is required to be 7 < Pt < 50 GeV. Events with at least two (three) jets with transverse 
momentum larger than 5 GeV are considered as dijet (trijet) events if the two leading jets of the 
measured observables have an invariant mass Mjj exceeding 16 GeV. For events that fulfill the 
dijet (trijet) criteria, the average transverse momentum is (Pt) = jj X]f Pp' *! with N = 2(3). 
All three double differential jet measurements are normalised to the inclusive NC DIS mea- 
surement as function of Q^. The advantage of the normalised jet cross sections compared to 
absolute jet cross sections [2] are reduced systematic uncertainties on the experimental as well 
as on the theoretical side. 

3 Detector correction using regularised unfolding 

Due to kinematic migrations because of resolution and other effects and due to the limited 
acceptance of the detector, the data have to be corrected. For this purpose a multidimensional 
regularised unfolding procedure, including all correlations, is applied. This procedure makes 
use of a migration matrix A that connects the particle level, represented by the vector x, with 
the detector level, represented by the vector y, such that the equation y — Ax holds. The 
particle level distribution x is determined using the TUnfold package [B]. The x^-function 

xH^) = {y-A-xf V-y' iy-A-x) + t\x - Xof{L^L)ix - Xo) (1) 

is minimised analytically as a function of x, where Vyy is the covariance matrix, and xq is 
the bias distribution. The Tikhonov regularisation parameter t protects the result from large 
fluctuations, and a regularisation condition for the matrix L the unit matrix is chosen. 

The unfolding is performed using a single matrix A with an overall 4x4 structure that allows 
to unfold four measurements (NC DIS, inclusive jets, dijet, trijet) all at once and further gives 
the possibility to determine the normalised cross sections taking all correlations into account 
(see sec. |4]). The four diagonal elements are submatrices that describe the migrations of just 
one single measurement. 

NC DIS submatrix: The migrations of the NC DIS measurements are described in two di- 
mensions, i.e. in the kinematic variables y and Q^. The NC DIS data are used to determine 
the normalised jet cross sections. 

Inclusive jets submatrix: The inclusive jet measurement is a measurement of jet multiplici- 
ties, where jets are defined by the jet algorithm independently on particle level and on detector 
level. This implies that jets have to be connected between both levels. A geometric jet matching 
method based on a closest pair algorithm with a distant measure R = ^A02 + Arj^ and a max- 
imum distance of P < 0.9 is applied. This ensures that no kinematical biases are introduced. 
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Matched jets are filled into the sub- 
matrix of the inclusive jets in a 
three-dimensional unfolding scheme 
in the variables Pt, and y. Jets 
that appear on generator level only 
are treated as inefficiencies. Jets 
that appear on detector level only 
are difficult to handle and cause 
a large error on the overall nor- 
malisation if not handled properly. 
Detector-level-only jets do not have 
particle level jet quantities, but still 
event observables on particle level 
are known. These jets are filled 
into the submatrix that connects 
the detector level inclusive jets to 
the particle level NC DIS measure- 
ment. The normalisation of the NC 
DIS measurement is preserved by 
adding negative weights to the ef- 
ficiency accordingly. 
Dijet and trijet submatrices: 
The unfolding of the dijet and tri- 
jet measurement is performed in a 
three dimensional unfolding scheme 
in the kinematic variables (Pt), 
and y. Additional bins describe 
migrations in and out of the dijet 
and trijet phase space. Migrations 
in the kinematic variables Mjj > 
16 GeV, P^*^^ > 5 GeV and require- 
ments on r]l^ are taken into ac- 
count. Similarly as for the inclu- 
sive jets, events that do not ful- 
fill the dijet (trijet) requirements 
on particle level but on detector 
level are estimated by the NC DIS 
events. 
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Figure 1: Normalised inclusive jet, dijet and trijet mea- 
surements in six bins of the virtuality of the exchanged 
boson and in bins of the transverse jet momentum Pt 
in the Breit frame or the average transverse momentum of 
the two (three) leading jets (Pt), respectively. The corre- 
lated statistical errors are shown by the inner error bars 
which are mostly smaller than the marker size. The outer 
error bar shows the total experimental uncertainty, incor- 
porating all systematic experimental uncertainties. The 
data are corrected for detector effects using a regularised 
unfolding technique. The measurements are compared to 
next-to-leading order QCD calculations using the CTIO 
PDF set. They are corrected for hadronisation effects. 
The theory predictions include uncertainties determined 
by scale variations of a factor 2 up and down. 



The covariance matrix on detector level Vyy is determined by data and contains the statis- 
tical uncertainty and all correlations between the four measurements. The unfolded covariance 
matrix Vxx on particle level is determined by error propagation through the unfolding pro- 
cess and holds information on the correlations resulting from detector effects as well as the 
propagated correlations between the single measurements. 
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4 Result 



4.1 Normalised multijet cross sections 

Each bin of the jet measurements is normahsed with the corresponding bin of the inclusive 
NC DIS measurement. The covariance matrix of all three normalised jet measurements is de- 
termined by a full error propagation using the covariance matrix V^x ■ Also the systematic errors 
are determined by a full error propagation. This procedure is equivalent to a direct measure of 
normalised jet cross sections. The normalised inclusive jet, dijet and trijet measurements are 
shown in Fig. [l] where they are compared to pQCD predictions which show a good agreement 
over the full phase space. 

4.2 Determination of the strong coupling constant 

The strong coupling is determined by performing a ^^-minimisation procedure to all three 
normalised jet measurements simultaneously with as{Mz) as a free parameter. The theory 
calculations are performed using the QCDNUM program [7] for the NC DIS cross sections and 
the NLOJET++ program [5112] interfaced to FastNLO [101 E] for a fast repeated calculation 
of the jet cross sections using the CTIO PDF set [T^. The x^-definition includes the full 
covariance matrix after the unfolding and takes into account systematic uncertainties using 
nuisance parameters. In order to consider only bins that show a fast convergence, bins with 
large fc-factors of fc > 1.3 are excluded from the fit, where k = ctnlo /'''lo ■ The resulting fit 
takes contributions from 42 out of 65 bins into account and shows a reasonable x^/ndf of 53/41. 
The resulting as{Mz) is determined to be 

as{Mz) = 0.1163 ± 0.0008 (exp.) ± 0.0011 (had.) ± 0.0014 (pdf) tomtt (theo.). 

This result is consistent with previous HI publications of normalised multi-jet cross sections 
[T3]. The uncertainties from hadronisation corrections (had.), the PDF uncertainty (pdf) and 
theoretical uncertainties from missing higher orders (theo.) are determined by repeating the fit 
with shifted theory cross sections. Each of the theoretical uncertainties is larger than the total 
experimental uncertainty. The result is consistent with the world average within the errors [14| . 

References 

[1] Hi Collaboration, preliminary result: Hlprelim-12-031 (2012) . 

[2] HI Collaboration, preliminary result: Hlprelim-11-032 (2011) . 

[3] R. Kogler. PhD thesis, Univcrsitat Hamburg, DESY-THESIS-201 1-003, MPP-2010-175, 2011. 

[4] S. D. Ellis and D. E. Soper. Phys. Rev. D 48 (1993) 3160. 

[5] M. Cacciari, G. P. Salam, and G. Soyez. Eur.Phys.J. C72 (2012) 1896, arXiv: 1111 .6097 [hep-ph]] 

[6] S. Schmitt. |arXiv : 1205 .6201 [physics ■ data-an]" 

[7] M. Botje. Comput.Phys.Commun. 182 (2011) , arXiv : 10 05 ■ 1481 [hep-ph]| 

[8] Z. Nagy and Z. Trocsanyi. Phys. Rev. D 59 (1999) 14020. 

[9] Z. Nagy and Z. Trocsanyi. Phys. Rev. Lett. 87 (2001) 82001. 

[10] M. Wobisch et al. arXiv: 1109. 1310 [hep-ph]' 

[11] D. Britzger et al. these proceedings (2012) . 

[12] H.-L. Lai et al. Phys.Rev. D82 (2010) 074024, arXiv: 1007.22 41 [hep-ph] [ 

[13] F. Aaron et al. Eur.Phys.J. C65 (20 10) ,,arXiv: 0904. 3870 [hep-ex] [ 

[14] S. Bethke. Eur. Phys. J. C 64 (2009) 689. 



4 



DIS 2012 



